NASA CR- 132625 


(NASA-CB-132625) COE CONTAMINATED N75-29562 

ATBO^THEBE: THE EAfcGER OF ClIHATE CHANGE, 

PHASES 1 ARE 2 Final Report (Crexel Oniv.) 

76 p HC J«.75 CSC I 13E Unclas 

G3/U5 31372 



EARTH ENERGY EXPERIMENT (E 3 ) PROJECT 

NASA CONTRACT NAS 1-11871 


drexel university 

PHILADELPHIA, PENNSYLVANIA 19104 




DEPARTMENT OF PHYSICS 
AND ATMOSPHERIC SOENCE 



NASA CR-132625 


OUR CONTAMINATED ATMOSPHERE - 
THE DANGER OF CLIMATE CHANGE 


By Alan J. Cimorelli and 
Frederick B. House 


Final Report 
Phases I and II 


Prepared under Contract No. NAS1-11871 by 
Departpert of Physics and Atmospheric Science 
Drexel University 
Philadelphia, Pennsylvania 19104 

for 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 



FOREWORD 


The present report is a Master's research paper prepared for the 
Environmental Science Department at Drexel University. This effort was 
partially supported by the Langley Research Center under Contract No. 
NAS1-11871, whose effort is directed towards an investigation of ESSA VII 
Satellite radiation data for use in long-term earth energy experiments. 

This report is one of three companion reports which together constitute 
the final report of phases I and II of subject contract. The other two 
reports are entitled; "An Investigation of ESSA VII Radiation Data for 
Use in Long-Term Earth Energy Experiments," published as NASA CR-3 32623, 
and "Steady-State Solution to the Conduction Problem of a Spherical Balloon 
Radiometer," published as NASA CR-132624. 



ABSTRACT 


Changes in average global and hemispheric surface temperature are 
related to climate change. The effects of increased concentrations of 
atmospheric particulate matter on average surface temperature and on the 
components of the earth's radiation budget are studied. To evaluate 
these effects. Sellers' (1973) model, which couples particulate loading 
to surface temperature and to changes in the earth's radiation budget 
was used. A determination of the feasibility of using satellites to 
monitor the effect of increased atmospheric particulate concentrations 
is performed. The principal findings of this study are: 

1) A change in man-made particulate loading of a factor of 4 is 
sufficient to initiate an ice-age. 

2) Variations in the global and hemispheric weighted averages 
of surface temperature, reflected radiant flux and emitted 
radiant flux are non-linear functions of particulate loading. 

3) Zonal variations in reflected and emitted radiant fluxes are 
linear functions of particulate loading in mid-latitudes and 
non-linear functions in higher latitudes. 

4) In order to detect changes in man-made particulate loading 
of a factor of 2 it was found necessary to sense changes in 
reflected and emitted radiant fluxes of 2.6 and 0.8 w/m^. 



respectively. This corresponds to changes in tlie equilibrium 
temperatui e of the monitoring instrument (black sphere) of 
-0.2°K (a sensitivity of one part in 1000) for night time 
measurements and -K).1°K (a sensitivity of 1 part in 3000) 
for day time measurements. 

It was found that a black satellite sphere proposed by House 
and Sweet meets the requirement of night time measurement 
sensitivity, but falls short of the required day time 
sensitivity. However, non-black, spherical radiometers 
whose external optical properties are sensitive to either 
the reflected radiant flux or the emitted radiant flux 
can meet the observational sensitivity requirements. 
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I. INTRODUCTION 


In recant years there has been much concern voiced as to the effect 
of man's activities on his health and his environment. Most of the 
effect, thus far, has been concentrated in those areas which directly 
affect his local environment; namely, the effect of increased ambient 
concentrations of SO 2 , N0 X , 002, 00, particulate matter, toxic trace 
metals, etc. The effect of such pollutants on man's health, and the 
local ecology has been researched to the point where criteria are 
available by which standards for allowable ambient concentrations 
can be set. 

Such standards are useful only if these concentrations can be 
effectively monitored so that a potential problem can be detected 
and corrected before it becomes a reality. For such local pollutant 
problems, effective ambient and source monitoring systems have been 
developed. Most techniques involve taking a small sample of air and 
chemically analysing it in order to determine relative pollutant con- 
centrations and then generalizing to a small local environment. 

Although these local pollutant problems are a threat to man’s well- 
being, they are reasonably understood and effectively monitored. 

Such environmental problems, at least for the time being, manifest 
themselves in enforcement. 

The area in which comparatively little concrete data exists for 
both determining environmental effects and developing reasonable 



monitoring techniques is man's effect on the global environment. Here 
there is a great deal of somewhat questionable evidence which suggests 
that due to man's activities, changes in the global environment may 
occur which will drastically threaten our life styles, if not our very 
existence. Specifically, has man inadvertently changed the global 
climate, or is he about to do so? 

It is necessary, according to Lorenz (1970) to first determine 
whether or not climate changes have occurred, and if so, their nature 
and extent. In order to deal with this problem in a complete way, it 
would be necessary to develop a means by which we could identify small 
trends in the characteristics in a widely random varying atmosphere. 
Although such a complete solution is not presently available, the 
evidence of climate change exists. History has recorded vegetation 
changes, in a particular locale, of a type which must have been caused 
by a dramatic change in rainfall or tempterature. For example, the 
great oil fields of Alaska, which were produced by vast amounts of 
organic matter, indicate that at the time of their inception, a climate 
much different from that of the present must have existed. Most 
striking is the advance and retreat of the prehistoric continental 
glaciers. These glacial epochs indicate large oscillations in the 
global climate, having periods in the order of millions or years. 

Aside from such long-term changes, there is flso evidence which 
indicates climate changes occurring in contemporary history. Mitchell 
(1971) among others, claims that during the ^ast century a systematic 
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fluctuation of global climate is revealed by meteorological data. He 
states that between 1880 and 1940 a net wanning of about 0.6°C occurred, 
and from 1940 to the present our globe experienced a net cooling of 
0.3°C. 

Once such changes have been established it is necessary that their 
causes and effects be identified. Towards this end a large number of 
variables which describe our climate must be identified and their 
relative importance assessed. Since this is such a large task, models 
which grossly parameterize these variables are developed, solved and 
perturbed. 

The effect occurring due to a climate change could be as mild as 
a slight inconvenience caused by hotter summers or as devastating as 
the annihilation of northern civilizations by the intrusion cf massive 
continental ice sheets. The effect is purely a matter of the magnitude 
of the climate change. 

The possibility exists, according to some authors, that our 
climate is intransitive. That is, small perturbations occurring in 
certain atmospheric variables may cause our climate to change to a 
new steady state climate. Quantities such as the amount of solar 
radiation received by the planet per unit area per unit time (solar 
irradiance), the total reflectivity of the earth-atmosphere system 
(albedo), and the concentrations of atmospheric constituents which 
affect the radiative transfer characteristics of the atmosphere 
(specifically} CO^ and particulate matter) are variables of this type. 



Immediately the question arises, is man in a position to affect these 
quantities and if so to what degree and with what resulting effect? 
Obviously man has no control over the incident solar irradiance. He 
can but monitor it and speculate on possihle climate change which would 
occur if it varied. Fortunately, this solar radiation averaged over a 
period of time is constant and there is little evidence which would 
suggest its changing. In recent years there has been much data collected 
which indicates that the global concentrations of both 00^ and parti- 
culate matter, through man's industrial activities, are on an increase. 

The effect of increasing global concentrations of CC >2 is to increase 
the average global surface temperature. This effect is known as the 
"Greenhouse Effect." It was first suggested by Tyndall in 1863, as 
reported in the SCEP (Study of Critical Environmental Problems) Work 
Group, that 

"... the hlanketing effect of increas^ 00 would cjuse 
climate changes through variation of the surface tempera- 
ture. Increased CO 2 , because of its strong absorptior. 

(and therefore emission) of infrared radiation at 12 to 
I8y, would reradiato surface and further inhibit the 
radiative cooling at the ground." 

As for increasing global concentrations of particulate matter, 
the opposite is true. Due to the backscattering characteristics of 
the particles an increase in the albedo occurs which produces a 
decrease in average global surface temperature. 

NASA Langley Research Center in the report "Remote Measurement 
of Poi i«+? on" states, 



"Particles affect the transfer of radiation through the 
atmosphere hy virtue of their ability to absorb, emit, and 
scatter radiation. In the visible portion of the spectrum, 
the scattering of radiation by particles can equal or exceed 
the scattering due to the normal gaseous constituents. In 
some cities the long-term average of solar radiation reac h i n g 
the ground is reduced by particles by more than l(#. n 

Since it had been thought that the two effects canceled one another, 

their use as criteria for developing standards has not been considered. 

However, it has since been found that the rate of temperature increase 

decreases with increasing CC>2 and increases with increasing particulates. 

Therefore, global particulate loading is of foremost concern. 

Sellers (1973) has developed a climate model which quantitatively 

relates particulate loading to surface temperature. He has stown that 

an increase in man-made global particulates by a factor of 4.0 will 

initiate an ice-age. In order that we safeguard ourselves and future 

generations from a self-imposed ice-age it is necessary that we 

effectively monitor globed concentrations of particulate matter. 

Using Sailers' data a determination of the feasibility of monitoring 

the effect of increases in global particulates through the use of existing 

satellite technology was performed. It was found that earth radiation 

balance measurements performed by a fleet of proposed satellites are 

accurate enough to implicitly detect, through earth radiation budget 

measurements, increases in global particulate concentrations. 
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II. BACKGROUND 


2.1 Earth's Radiation Budget 

Consider the earth and atmosphere as an entire system, then at the 
interface between the earth's atmosphere and space, incoming and out- 
going radiation fluxes are the means by which energy is exchanged 
between the sun and the earth - atmosphere system. The earth receives 
short wave radiation (.2y to 3.5y) from the sun (insolation). A portion 
of this radiation is reflected back into space and therefore contributes 
nothing in the way of available heat energy to the earth-atmosphere 
system. The remainder of the radiation is absorbed by the earth and 
atmosphere. This absorbed radiation acts to heat the absorber which 
in turn reradiates according to Stefan-Boltzmann ' s Law and at a 
wavelength consistent with Plank's Law. The earth loses energy 
(emitted radiation) to space by thermal emittance of longwave radia- 
tion at wavelengths greater than 4«0y. This rate of heat loss must 
exactly balance over a long period of time the amount of heat received 
from the sun, in order to prevent the earth from beco ming continuously 
hotter or continuously colder, as the case may be. This energy 
exchange with space is the radiation balance or net radiation at the 
top of the atmosphere. Linked to the radiation balance, through 
thermal exchange mechanisms, is the surface temperature of the earth. 

It is this parameter which is quantitatively synonymous with climate. 
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One other pertinent quantity needs to be defined; namely, the 
albedo . The albedo is the percentage of the insolation which is 
reflected back into space, i.e. the ratio of the reflected radiation 
to the insolation. 

Figure L shows the mean annual heat budget (radiation balance) of 
the northern hemisphere. Here is seen the way in which radiation dis- 
tributes itself within the earth-atmosphere system. The actual quantities 
appearing in the figure are Representative but somewhat out-dated. For 
example, recent satellite measurements indicate that the albedo is 29 
or 30Z. 

The earth's net radiation budget for any space and time scale is 
given by: 

Q - H - (W + W ) 
n s re 

or, Qn * H s (1“A) - W g 

where, H = H-. cos 0 

s * 

Figure 2 describes this situation pictorially. 

Satellite experiments are designed to effectively measure emitted 
radiant flux and albedo on a predetermined space and time scale. In 
order to generate the global net radiation, the reflected radiant flux 
is integrated over the entire earth. 

2.2 Global Climate Models - Developmental Summary 

In order to relate the quantities found in the radiation budget 
to surface temperature we must specify all those characteristics and 


- 7 - 




O 

i-l 

M 

a) d 
*C O 
0.1-1 



O 4J 

+ 

m 

4J <d 

m *h 


¥ 

<u 

o 

td 

O* 

00 

4J cd 
in m 

o 


O 

•H 

M 


4J 

o> d 

>n 


O 

• 

d 

O. 1-f 


0 

09 i-> 


•H 

o cd 

+ 

iJ 

O. tH 


cd 

o ♦d 
M d 


'g 

3 

H M 



d 

o 

t3 i-l 



d 4J 

*4 


3 cd 



O i-» 

+ 


M -d 
oo cd 

W 



CO 

<D 3 
4J O 
o f-H 
<U Cj CM 

*"H + 

m -4 a 
a) o 
M U 


T3 

0) 

V* 

0) 


c 

o 

a 



3H3HdS 

-oivaxs 


aaaHdsodOHi 


QNnOHO 


8 


Figure 1: The mean annual heat budget of the Northern Hemisphere, after London. 
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mechanisms which affect heat transfer within the earth-atmosphere system, 
i.e. , atmospheric motions, ground conditions, radiative and connective 
processes within the atmosphere and atmospheric constituent concentra- 
tions. To accomplish this, global climate models are developed. Using 
the data produced from earth radiation budget measurements as input for 
a global climate model, surface temperature fields can be generalized. 
The worth of these models is seen in their ability to reproduce our 
existing climate using available earth radiation budget data, whereas 
their merit is in their ability to predict climate change. 

In general, global climate models preceed as follows: In order to 

deal with climate on a mathematical basis, the climate must first be 
defined as a collection of long-term statistical properties of the 
atmosphere. Usually five such properties are used: the three dimen- 

sional fields of temperature, pressure, density, wind velocity and water 
in its various forms. Once these properties have been identified the 
physical laws which govern their variations must be written as mathe- 
matical equations. A solution of this set of equations defines the 
state of the atmosphere at a given instant of time. But, before a 
solution can be generated the state of the environment must be 
identified. This can be done in one of two ways. It can be assumed 
that a knowledge of how the environment will behave is known, in which 
case environmental factors enter the model as constants. Or we 
could introduce those characteristics which describe the state of the 
environment as variables in terms of them. For example, allow 
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particulate concentrations to enter the model in an equation which 
relates available solar radiation to environmental factors. 

To illustrate this consider the following example. In any globed 
climate model a fundamental quantity which must be determined in order 
that it be used as input to the model is the amount of solar radiation 
absorbed by the earth-atmosphere system. This quantity is found by 
differencing the solar radiant Hix incident upon the atmosphere, which 
is known to be approximately constant (solar constant), and the solar 
radiant flux reflected by the earth-atmosphere system. Moreover, the 
reflected radiant flux depends strongly upon environmental conditions 
such as cloud cover, type of underlying surface and concentration of 
atmospheric constituents. Specifically, it is well known that the 
reflected radiant flux is sensitive to atmospheric particulate con- 
centrations; therefore, our climate model must account for such 
concentrations. This can be accomplished in one of the two ways 
previously stated, where in this case the environmental factor is 
atmospheric particulate concentrations and the atmospheric variable 
is the reflected radiant flux. Once the equations are written and 
environmental factors considered a solution proceeds as follows: 

1) In the equations all partial derivatives are replaced by 
finite differences. 

2) The initial and boundary conditions are chosen. 

3) Special and temporal increments are chosen. 

4) The equations are then solved in a step-wise manner with 
the aid of a digital computer. 
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Since there is such a large number of variables to be considered, 
those models which are capable of being completely solved are necessarily 
quite simple. 

2.3 Global Climate Models - Futuristic Pro .lections 

Once it has been established that a climate model accounts rea- 
sonably well for the present climate it is then used as a speculative 
tool. This is accomplished by purturbing certain parameters used in 
the model and examining the resulting climate change. Although such 
models differ in their treatment, their results are quite similar. 

Budyko's model (1970) indicates that a reduction of 1.6$ of the 
solar radiation will produce an ice cover on the globe reaching a mean 
latitude of 50°. Sellers (1973) found a similar result using his 
model . He determined that a decrease of 2$ in solar radiation wo . 
produce ice caps extending to 50°. Most frightening of all is the woi 
done by Faegre (1972) : he reports that there are five steady state 

solutions (global temperature distributions) consistent with his 
model. One of the solutions corresponds to our present climate while 
the others represent much different climates. In perturbing each of 
the solutions he found that our present climate may be unstable and 
that a decrease of 2$ in radiation results in a climate corresponding 
to an ice-covered earth. Although this result is consistent with the 
results of Budyko and Sellers, Faegre 's model is dissimiliar in that 
it indicates that if the earth is perturbed into an ice-covered climate 
it will not return to our present climate should the radiation return 
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to its present value. Even though these models contain a large number 
of assumptions which grossly simplify the actual systems and for this 
reason some scientists refute them, they are the best working tool 
available in the area. Since there is such close agreement among 
them, it seems quite necessary to determine which atmospheric consti- 
tuents and or pollutants, controllable by man, are able to affect such 
climate change. Having identified these constituent^ how can we 
monitor and control their numbers? 


2.4 Effect of Particulates and Carbon Dioxide on Climate 


It has been well established, in the literature for a number of 
years, that the two pollutants which possess the greatest possibility 
of drastically effecting our climate are 00^ and particulate matter. 
Up until recently, it was thought that the cooling effect produced by 
increased concentrations of particulate matter would be canceled by 
the warming effect of CC^. But it was found by Rasool and Schneider 
in 1971 that, 

”... although the addition of CXD^ i n the atmosphere 
does increase the surface temperature, the rate of 
temperature increase diminishes with increasing CO5 in 
the atmosphere. For aerosols, however, the net effect 
of increase in density is to reduce the surface tempera- 
ture of the Barth. Because of the exponential dependence 
of the backscattering the rate of temperature decrease is 
augmented with increasing aerosol content." 

Their results showed that an increase by a factor of eight of global 

CO2 concentrations will produce an increase in surface temperature 

of less than 2°C, whereas increasing particulate concentrations by a 
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factor of four could decrease the mean surface temperature as much as 
3.5°C. Therefore, as the global concentrations of both CO2 and parti- 
culates are increased, it is the radiative shielding effect of the 
particulates which plays the dominate role.- 

Although such an increase in particulate pollution seems quite 
unlikely, a number of studies indicate that such an occurrence is 
possible if we continue polluting at our present rate. In a pubication 
prepared by NASA Langley Research Center (1971) continued increasing 
particulate concentrations are reported, and I quote, "There is little 
doubt that over many populated areas, atmospheric turbidity is 
increasing. Solar radiation measurements in Japan and in the European 
part of the USSR show this increase clearly.” MacCormick and Ludwig, 
1967, presented data which indicated that global turbidity has 
increased approximately 1.0 to 2056 per decade for the general period 
of 1910 to I960. Their data comes from measurements taken from Davos, 
Switzerland and Washington, D.C. Also, Schaefer's, 1968, measurements 
taken in the United States, have shown that atmospheric particulate 
concentrations have increased by at least an order of magnitude in 
the last ten years. Finally, Ludwig and Morgan, 1970, predict that 
man's potential to pollute will increase six-to-eight fold in the 
next fifty years. All of these studies point to the possibility, in 
the not so distant future, of man polluting himself into an ice age. 

Sellers, 1973, developed a new global climate model which 
accounts for significant atmospheric variables. He found, upon 
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analyzing the sensitivity of his model to increased concentrations of 
OO 2 and particulates, a similar result to that of Rasool and Schneider. 
He doubled and halved the initial CCU content of the atmosphere with no 
appreciable effect on surface temperature. However, upon analyzing the 
model 1 s sensitivity to changing concentrations of particulate matter he 
found surface temperature to drop even more drastically than that pre- 
dicted by Rasool and Schneider. His model predicts a drop of 0.5°C in 
global surface temperature by doubling man-made global atmospheric 
particulate concentrations. And if we increased the amount of man- 
made particulates by a factor of 3 global surface temperature would 
decrease 1.3°C. His findings show that the rate of decrease of global 
surface temperature increases as more particulates are added to the 
atmosphere. He increased the amount of man-made global particulates 
as much as six fold which resulted in a global surface temperature 
drop of 12.5°C. If man is in a position to affect such change then 
it is certainly the tine to reevaluate our priorities. 

2.5 Monitorin g Partic ulate T.oading 

Unlike pollution problems on a local level, the problem of 
increasing particulate concentrations on a global basis is quite 
vague in comparison. 

The studies and models cited above, by their very nature, do 
not afford working criteria by which reasonable standards could be 
set and control implemented. But, since they do indicate the 



possibility of such drastic global consequences the implementation of 
some form of early warning monitoring system is essential. 

The only methods available, by which a determination of the ambient 
concentrations of particulates can be made, are good only for small 
local areas. Such data is useless for monitoring climate change due 
to increased particulate concentrations since average global concentra- 
tions are needed. Therefore, rather than monitoring these concentrations 
directly, we could monitor some quantity whose magnitude is quite sensi- 
tive to changes in particulate concentrations. One way to accomplish 
this would be to monitor changes in average global surface temperature 
since according to Sellers (1973) changes in surface temperature are 
certainly sensitive enough ( see figure 3 and 4) . However, the problem 
with using such a quantity is that our surface sensing devices do not 
have global extent. And even if an accurate measure could be made by 
the time we sensed the change it would be unlikely that anything could 
be done about it. Therefore, a parameter is needed which is first, 
sensitive enough to changes in atmospheric particulate concentrations, 
such that its resulting change is within the limits of its measure- 
ment, and secondly, can be measured on a global basis. Earth energy 
budget measurements were found to be of this type. 
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III. APPROACH 

Sellers 0.973 ) reported the variation in average global surface 
temperature as a function of particulate optical depth. It was neces- 
sary tc relate these data to what occurs at the top of the atmosphere. 
That is, those changes in the components of the global radiation budget 
which must have occurred to produce the reported changes in the average 
global surface temperature. In order to accomplish this a request was 
made to Dr. Sellers for the data which were used as input to his model 
to generate the surface tenperature variations. Whereupon data were 
sent concerning zonal variations in net radiation, albedo as a function 
of molecular plus particulate optical depth for January and July, and 
the variation in solar irradiation as a function of latitude for the 
twelve months. 

Using these data it was possible to determine the sensitivity of 
the short and long wave radiant exstance to changes in global parti- 
culate concentrations. Given the net radiation, albedo, solar irra- 
diance for a specific latitude zone, optical depth and time of year, 
the corresponding emitted and reflected radiant fluxes were calculated 
as follows: 

w e i3k = V k d-A ijk ) - On 1 ” ' (!) 

and 

W r iJk = H 0 Jk A^ k (2) 
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The sensitivity of reflected and emitted radiant fluxes to cha’igeo 
in the global concentrations of particulate matter was compared to a 
satellite's sensitivity to changes in the reflected and emitted radiant 
fluxes. This was done in order to determine the feasibility of using 
satellites to monitor the effect of increased global concentrations cf 
particulates. 

Again, using SelLers' data it was possible to compute the global and 
hemispherical area-veLghted averages for solar reflected radiation and 

terrestial emitted radiation for January and July for varing particulate 
concentrations. The averages were calculated as follows: 



h=o 

In order to determine whether or not these changes in the earth 
radiation budget components with particulate loading can be sensed by 
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the satellite monitoring instrument, it was necessary to compute the 
changes In equilibrium temperature of the Instrument with particulate 
loading. It is the equilibrium temperature of the instrument which is 
measured and from which earth radiation budget components are deciphered. 
These changes in equilibrium temperature were then compared with last. 
went sensitivity in order to determine the feasibility of using such a 
monitoring system. 

The instrument considered is a metallic spherical detector optically 
black. The external radiation absorbed by the sphere (the addition of 
solar incident, solar reflected and terrestrial emitted) produces a 
proportional Internal irradiance that will- be detected, as equilbrlun 
temperature, by internally facing r a di o m e ters. 

The equations used to c alc u late the equilibrium temperature of 
the sphere are as follows: 

When the sphere is illuminated by the sun (daytime) 

* "e + (5) 

When the sphere is in the shadow of the earth (nighttime) 

T " C (P a-e V 1 a J 1/4 < 6 > 



.2705 for a sphere at 800 km assuming 
isotropic conditions. 
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IV. RESULTS 


Through the use of the climate model developed by Sellers (1973) 
the feasibility of monitoring climate change, i.e. changes in surface 
temperature, due to increases in gobal atmospheric particulate concen- 
trations is reported. 

It was stated previously that a change in average surface tempara- 
tuer of 3°C is sufficient to initiate an ice age. If the monitoring of 
climate change is to be a viable tool through which global atmospheric 
concentrations can be controlled it would be necessary for the monitoring 
instrument to be able to sense changes in average surface temperature 
of much less than 3°C. Changes in the order of one-half degree would 
be meaningful. 

Using the data reported in Appendix A global and hemispheric 
weighted averages and zonal variations in surface temperature, short- 
wave reflected radiant flux and long-wave emitted radiant flux for 
January and July as functions of molecular plus particulate optical 
depth are generated. The results of these computations are found in 
Appendicies B and C. These results are also plotted in Figures 3 
through 14* 

On the graphs, an optical depth of 0.300 corresponds to the present 
optical thicknesb of the clear atmosphere. This consists of the sum 
of a molecular optical depth of 0.145 and a particulate optical depth 



of 0.155. Of the particulate optical depth it is thought that man con- 
tributes 0.023, or about 15%. The optical depth was increased by multiples 
of 0.023 such that upper limit of particulate loading, i.e. 0.415, 
appearing on the graphs corresponds t* ■ an increase in man-made parti- 
culate optical depth by a factor of s x. 

4.1 Zonal Variations in Earth Radiation Budget Components with Particulate 

Loading 

Zonal variations in reflected and emitted radiant fluxes and albedo 
with particulate loading for January and July are computed. The results 
of these computations are found in Appendix B. These results are also 
plotted in Figures 3 through 8. 

Figures 3 through 8 show the variations reflected and emitted 
radiant fluxes and albedo with particulate loading for January and July. 

The variations are calculated for each 10 degree latitude zone. In 
the figures the latitude zones of 60-70 north, 30-40 north, 0-10 north, 
30-40 south and 60-70 south are reported. 

A general characteristic of the graphs is that in the low and 
middle latitudes the dependence of th< reflected radiation stream and 
albedo on particulate loading is essentially linear whereas the depen- 
dence of the emitted radiant flux is i.on-linearal 1 y decreasing. In 
the higher latitudes the dependence of the reflected and emitted 
radiant fluxes are found to be non-linearally increasing and decreasing, 
respectively. These results can be accounted fcr in the fact that for 
high latitudes the incident solar radj ation incounters a larger atmos- 
pheric path than in the low and mid latitudes. Therefore, the radia- 
tion stream comes in contact with a larger number of particles which in 
fact increases the effective optical depth at those latitudes. 



Particulate Loading (optical depth) 

Figure 3 . Zonal Reflected Radiant Flux for January vs. Particulate Lo ading 
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P Articulate loading (optical deplih) 

Figure 4* Zonal Reflected Radiant Flux for July vs. Particulate Loading 
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Particulate Loading (optical depth) 


Figure 5* Zonal Emitted Radiant Flux for January vs. Particulate Loading 
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Particu ate L odi (optic-r 5 ! depth) 


Figure 8. Zonal Albedc for July vs. Particulate Loading 
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4.2 Variations In Global and H wnd anhar jc Weighted Averag es nf Ewpth 

Radlat 4 on Budget Componatits with Particulate Loading 

Variations In global and hemispheric weighted averages for surface 
temperature and reflected and emitted radiant fluxes with particulate 
loading for January and July are computed. The resalts of these com- 
putations are found in Appendix C. These results are also plotted in 
Figures 9 through 14* 

Figures 9 and 10 show the variation in average surface temperature 
with particulate loading. It is found that their functional relation- 
ship is non-linear ally decreasing. This agrees well with Rasool and 
Schneider in that the rate of temperature decrease is augmented with 
increasing aerosol content. This non-linear dependence is due to the 
backscattering characteristics of the particles. It can be seen by 
examining the graphs that in order for t'le average global surfa ;e 
temperature to decrease by 3°C man-made particulate loading must 
increase by a factor of 4. This increase If seen as an incroare in 
tc *l atmospheric optical depth from 1 .300 to 0.369* Ap previously 
stated our interest is in the ability to ietect a decrease in surface of 
one-half degree. Examining the graphs re/eals that for such a decrease 
to occur man-ma.de particulate loading must increase by a factor of 2, 
corresponding to an increase in total atmospheric optical depth from 
0.300 to 0.323. 

Figures 11 and 12 show the variations in the global and hemis- 
pheric weighted averages of short-wave reflected radiant flux with 
particulate loading for January and July. It is seen that the 
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Particulate Loadi i.r (opticrl depth) 


Figure 10. Global and Hemispheric Average Surface Temperatures 
for Judy v.s. Particulate Loading 
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Particulate Loading (optical depth) 


Figure 11. Globed, and Hemispheric Average Reflected Radiant flux 
for January v.s. Particulate Loading 
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functional relationship between average reflected radiation and parti- 
culate loading is non-1 inear ally increasing. Since the incident 
short-wave reflected radiant flux which is reflected back into space is 
not available for absorption by the earth's surface there exists an 
inverse relationship between surface temperature and reflected 
radiation. Is reflected radiant flux changes surface temperature 
changes in the opposite sense. Therefore, since the functional 
relationship between surface temperature and particulate loading 
is non-] inearal ly decreasing the result cited here is reasonable. 

Upon m—in-ing tables C5 and C6 it is seen that in order to sense 
a change in particulate optical depth of 0.023, our monitoring instru- 
ment must be able, on the average, to sense a change in reflected 
radiant flux of approximately 1.3 watts/meter . 

Figures 13 and 14 show the variations in the global and hemispheric 
weighted averages of long-wave emitted radiant flux with particulate 
loading for January and July. Similarily, it is semi that the 
functional relationship between average emitted radiant flux and 
particulate loading is non-1 in earally decreasing. This is due to the 
direct relationship which exists between emitted radiation on surface 
temperature. The amount of emitted radiant flux is determined in part 
by the tenperature of the emitting surface. This relationship is known 
as St ef an-Bol t zmann 1 s Law, i.e. emittance is proportional to the forth 
power of the tenperature of the emitting surface. 

Again we wish to be able to detect changes in average emitted 
radiant flux which corresponds to a change in particulate optical 
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.tart ictl ate r ,oadi/;'' (optical depth) 


Figure 13* Global and Hemispheric Average Quitted Radiant Flux 
for January v.s. Particulate Loading 
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Figure 14. 



Particulate Laodinc (optical iep ;h) 


Global and Hemispheric Average Quitted Radiant Flux 
for July v.s. Particulate Loading 
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fro® 0.300 to 0.323* Engining Tables C3 and C4 reveals that in order 
to sense this increase in optical depth our monitoring instrument must 
detect a change in emitted radiant flux of approximately 0.8 watts/ 


4.3 The Feasibility of Monitoring the Effect of Increased Global 

Particulate Concentrations 

A determination of the feasibility of detecting changes in the 
components of the earth's radiation budget which correspond to a 
particulate loading increase of 0.023 is performed. In order to 
determine the feasibility of detecting such changes in the earth's 
radiation budget it is necessary to compute the changes in equilibrium 
temperature of the instrument with particulate loading (see approach) . 
The instrument considered is a spherical detector having an eoissivity 
equal to unity (black sphere). The results of these computations are 
given in the tables below. The global weighted averages for W Q and W r 
were used in equations (5) and (6) . 


j?a 


articulate] 

fading 
Bg * — y> , „„ 


.300 




A T 


ZgMlggflM 


1 .346 

AT 

11 ,?oo*3j 



.360 I AT 


TABLE 1: Change in Equilibrium Temperature of a Black Sphere as a 

Function of Particulate Loading for January. 



TABLE 2: The Zonal Variation of Net-Radiation with Particulate 
Loading for January 
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It is seen that day time measurements are much less sensitive to 
changes in particulate loading than are night time measurements. This 
is due to the fact that as W r increases with increased particulate 
loading V Q decreases ah amount such that the resultant change to be 
sensed by the instrument is nearly canceled. It is seen that a 
sensitivity of one part in 3000 is needed in order to sense a change 
in optical depth of 0.023. Therefore, by using a black sphere as the 
monitoring instrument all daytime measurements may be in the noise 
of the instrument. If one wishes to use daytime measurements a 
spherical satellite is needed whose optical characteristics are more 
sensitive to either V e or V r radiant fluxes. For example, an aluminum 
sphere is 3 to 4 times more sensitive to W r than is the black 
sphere. Similarly, white optical characteristics are such that 
they reflect in part and absorbs W g in a manner similar to 

a black sphere. 

Examining tables 1 and 2 (black sphere measurements) reveals 
that in order to be able to sense a change in average surface tem- 
perature of one-half degree or an increase in global atmospheric 
particulate concentrations of 0.023 a satellite monitoring system is 
needed which can sense a change in its equilibrium temperature of 
one part per thousand. Such measuring percision can be accomplished 
at night since there is no reflected radiant flux at this time and 
therefore changes occurring in the emitted radiation can be sensed by 
the instrument. 
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House and Sweet have proposed to _AAFE a long term zonal energy 
budget experiment employing a satellite network which is ideal for the 
particulate monitoring system described. In the proposal they report 
an ability to sense changes in the order of one part per thousand. 

Such precision would be adequate to monitor the effect on climate of 
increases in global particulates. 
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V. CONCLUSIONS 

Using data, which were generated from a global climate model 
denr&ped by Sellers, it was possible to determine the feasibility of 
monitoring the effects of increased global particulate concentrations 
through the use of existing satellite technology. It was determined 
that earth radiation budget measurements performed at night by a fleet 
of proposed Hack spherical satellites (House and Sweet) are accurate 
enough to sense a change in global particulate concentrations corres- 
ponding to a decrease in average global surface temperature of 0.5°C. 

Daytime measurements were found to be of use only if satellites, 
whose surfaces opticly discriminate between either W Q or W r , are used. 
It was pointed out that an aluminum or white satellite sphere could be 
used for the purpose of daytime measurements. It was also found that 
variations in the global and hemispheric weighted averages of surface 
temperature, reflected radiant flux and emitted radiant flux are 
non-linear functions of particulate loading. 

Since Sellers 1 model was used as the basis of the sensitivity 
calculations, the results can only be as good as the assumptions used 
in the model. According to Sellers, "Although the model is reasonably 
complete, in that most of the important atmosphere-land-ocean inter- 
actions are included, there is much room .. ir improvement of the various 
parameterizations." Therefore, these results should only be used to 
indicate the possibility of using such a technique for global 
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particulate concentration measurements. 

As a follow-up and or check of the validity of these results one 
could, with a complete record of spacial and temporal variation in 
albedo, attempt to correlate incidences of high particulate loading 
(e.g. volcanic erruptions with albedo change.) 
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APPENDIX A 

Zonal Variations in Solar Constant w:.th 
Month of Year and Zonal Variation In 
Net Radiation with Particu i«tft T-oarHnfl 

In this section is given, in tabular form, the data which was 
received from Dr. Sellers upon written request. These data include: 

1) the variation in solar constant with latitude for the twelve months 
of the year; (one month = 30.4375 days) j 2) the variation in net 

i 

radiation average surface temperature and albedo with particulate 
loading for each ten degree latitude zone for January and July. 

These data were provided in units of kly/So & kly respectively. 

It was found necessary, for reasons of comparison, to convert from 
kly to watts/m. Therefore, the data reported in this Appendix are 
in units of w/m . 
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TABLE A2: The Zonal Variation of Net-Radiation with Particulate 

Loading for January (Qq in watts/m^) 


latitude 

Sone 

Particulate Loading I 

.300 

.323 

.346 

.369 

^392 


?0-80N 

-163.0 

-160..6 

-158.7 

-154.2 

-147.1 

-135.1 

JO-70N 

-169.9 

-167.8 

-165.1 

-161.2 

154.4 

-151.5 

70-6GN 

-177.3 

-175.7 

-172.5 

-168.6 

-160.1 

-147.6 

dO-5QN 

-160.1 

-160.1 

-160.3 

-160.3 

-156.2 

-146.O 


-117.6 

-117.3 

-118.3 

-120.9 

-124.6 

-133.2 

40-30N 

-82.00 

-81.34 

-80.23 

-77.04 

-71.63 

-77.04 

3C-20N 

-44.41 

-44.10 

-43.46 

-40.75 

-35.97 

- 25.47 

20-10N 

-0.478 

-0.478 

0.000 

2.228 

6.367 

15.60 


36.93 

36.93 

3^.25 

39.32 

43.30 

52.05 

0-10S 

64.63 

64.47 

64.63 

66.70 

70.67 

79.27 

10-205 

87.07 

86.59 

86.75 

88.66 

92.32 

100.9 

20-30S 

100.6 

100.1 

100.1 

101.9 

105.5 

114.3 

30-40S 

114.3 

113.7 

113.7 

315.6 

119.4 

128.5 

40-50S 

114.4 

113.8 

114.0 

116.2 

120.3 

81.34 

50-60S 

98.85 

96.78 

91.05 

74.49 

42.34 

-15.60 

60-70S 

24.83 

23.08 

16.71 

0.159 

-25.79 

-49.98 

70-80S 

-49.50 

-49.19 

-51.26 

-59.37 

-54.60 

-45.52 

80-90S 

-54.76 

-54.12 

-52.53 

-49.34 

-44.89 

-36.29 








TABLE A3: The Zonal Variation of Net-Radiation with Particulate 

Loading for July ((^ in vatts/m z ) 


Latitude 

Zone 

Particulate Loading 

.300 

.323 

.346 

.369 

.392 

.415 

30-8QN 

-20.37 

-30.08 

-47.12 

-50.62 

-43*14 

-31.34 

3O-70N 

7.004 

-2.228 

-14.64 

-43.30 

- 57.46 

-44.25 

70-6QN 

29.75 

19.10 

1.592 

-37.25 

-58.90 

-56.35 

SO-50N 

76.09 

76.72 

78.31 

58.90 

5.730 

-51.10 

50-40N 

83.89 

84.05 

85.16 

88.98 

95.51 

19.58 

40-30N 

89.30 

89.14 

89.62 

92.64 

98.05 

110.2 

30-2CH 

73.22 

73.06 

73.38 

75.93 

80.86 

91.37 

20-10N 

58.10 

57.94 

58.42 

60.65 

65.10 

74.34 


43.30 

43.14 

43.46 

45.68 

49.82 

58.58 

3-10S 

17.67 

17.35 

17.67 

19.74 

23.88 

32.63 

L0-20S 

-15.92 

-16.24 

-15.92 

-13.85 

-0.710 

-0.637 

20-30S 

-58.74 

-58.74 

-58.26 

-55.87 

-51.57 

-41.86 

30-40S 

-95.51 

-95.35 

-94.71 

-92.00 

-87.39 

-77.20 

40-50S 

-13.18 

-131.6 

-130.7 

-12.75 

-125.1 

-131.5 

50-60S 

-162.4 

-162.4 

- 162.0 

-16.11 

-J. 6 O .4 

-153.6 

60-70S 

-178.8 

-178.4 

-177.2 

-17.30 

-166.7 

-154.4 

70-80S 

-171.3 

-171.0 

-169.5 

-16.52 

-159.0 

-147.4 

B0-90S 

-156.0 

-155.7 

-154.6 

-15.06 

-144.5 

-133.5 








TABLE A4> The Zonal Variation of Average Surface 

Temperature with Particulate Loading for 
January (Degrees K) 


Latitude 

Zone 

Particulate Loading 

.300 

.323 

.346 

•369 

.392 


90-8QN 

245.1 

243.6 

241.6 

238.7 

233.7 

224.9 

8O-70H 

251.1 

249.7 

247.8 

245.1 

239.8 

230.9 

70-6QH 

259.8 

258.7 

257.1 

253.7 

247.8 

238.5 

6O-50H 

268.9 

268.1 

266.9 

264.0 

259.1 

248.7 

5O-40M 

277.9 

277.3 

276.3 

274.0 

270.3 

262.0 

40-30H 

286.4 

285.8 

285.0 

283.1 

279.9 

273.3 

30-3CK 

292.2 

291.7 

290.9 

289.3 

286.4 

280.6 

20-1QH 

296.0 

295.5 

294.8 

293.2 

290.5 

285.2 

10-OH 

298.1 

297.6 

296.9 

295.4 

292.7 

287.5 

0-10S 

299.0 

298.5 

297.8 

296.2 

293.5 

288.3 

10-20S 

299.4 

299.0 

298.3 

296.8 

294.3 

289.2 

20-30S 

297.9 

297.6 

297.0 

295.4 

292.9 

287.8 

30-40S 

293.2 

292.9 

292.2 

290.6 

288.0 

282.5 

40-50S 

286.4 

286.2 

285.5 

283.7 

280.9 

274.7 

50-6QS 

278.9 

278.7 

277.9 

276.0 

273.0 

266.8 

60-70S 

271.8 

271.6 

270.9 

296.1 

266.3 

260.4 

70-80S 

266.7 

266.6 

265.8 

263.9 

261.2 

255.6 

60-90S 

262.9 

262.7 

262.1 

260.3 

257.7 

252.2 
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TABLE A5: The Zonal Variation of Average Surface Temperature 

with Particulate Loading for July (Degrees K) 


Latitude 

Particulate Loading 

Zone 

.300 

.323 

.346 

.369 

-•392 

.415 

90-808 

268.0 

267.3 

266.0 

262.8 

257.7 

249.3 

80-708 

272.9 

271.9 

270.4 

267.0 

261.8 

253.0 

70-608 

280.2 

279.1 

277.3 

273.1 

267.4 

258.0 

60-508 

286.5 

285.7 

284.4 

280.9 

275.3 

265.0 

50-408 

292.3 

291.6 

290.5 

288.0 

283.6 

274.7 

40-308 

296.7 

296.2 

295.3 

293.2 

289.8 

282.9 

30-208 

299.4 

298.8 

298.1 

296.2 

293.1 

287.1 

20-108 

299.3 

298.9 

298.2 

296.5 

293.6 

288.1 

LO-OB 

298.8 

298.4 

297.7 

296.0 

293.3 

288.1 

0-10S 

298.1 

297.7 

297.1 

295.5 

292.9 

287.7 


296.6 

296.2 

295.7 

294.1 

291.6 

286.4 

20-30S 

293.2 

292.9 

292.3 

290.8 

288.2 

282.8 

30-40S 

287.6 

287.4 

286.8 

285.1 

282.3 

276.3 

40-50S 

279.7 

279.6 

278.9 

276.9 

273.6 

265.3 

50-60S 

270.6 

270.4 

269.6 

266.9 

262.6 

253.0 

60-70S 

260.8 

260.6 

259.7 

256.5 

252.1 

243.3 

70-80S 

252.4 

252.2 

251.3 

248.5 

244.1 

235.7 

80-90S 

24^*4 

246.2 

245.3 

242.5 

238.2 

230.0 






TABLE A6: The Zonal Variation of Albedo with 

Particulate Loading for January (A injt) 



Particulate Loading 


.300 

.323 

.346 

.369 

.392 

.415 

90-80N 

74.5 

74.3 

74.2 

74.0 

73.8 

73.6 

80-70H 

73.0 

72.9 

72.8 

72.7 

72.6 

72.5 

70-6GH 

67.2 

67.7 

68.6 

71.1 

72.9 

72.8 

60-50B 

48.7 

50.6 

53.6 

60.0 

65.4 

72.8 

50-4QH 

32.8 

33.4 

35.4 

39.9 

47.0 

62.9 

40-30H 

29.4 

29.6 

29.8 

30.0 

30.2 

38.5 

30-20H 

27.9 

28.1 

28.4 

28.6 

28.8 

29.1 

20-10H 

25.9 

26.1 

26.4 

26.7 

26.9 

27.2 

imp 

26.0 

26.2 

26.5 

26.7 

27.0 

27.2 

p-ios 

26.7 

26.9 

27.2 

27.4 

27.6 

27.9 


26.5 

26.8 

27.0 

27.3 

27.5 

27.7 

20-30S 

26.9 

27.2 

27.4 

27.7 

27.9 

28.1 

30-40S 

27.2 

27.5 

27.7 

28.0 

28.2 

28.4 

40-50S 

30.3 

30.5 

30.7 

30.9 

31.1 

41.2 

50-60S 

34.7 

35.2 

36.6 

40.8 

47.9 

62.9 

60-70S 

51.9 

52.3 

53.9 

58.1 

64.6 

71.6 

70-80S 

70.1 

70.1 

70.8 

73.1 

73.0 

72.9 

80-90S 

74.7 

74.6 

74.5 

74*4 

74.2 

74.1 







TABLE A7: The Zonal Variation of Albedo as a Function 

of Particulate Loading for July (A in %) 


Latitude 

Particulate Loading 

Zone 

«30Q 

.323 

.346 

.369 

.392 


90-808 

63.6 

65.9 

69.9 

71.6 

71.5 

71.5 

80-701 

54.3 

56.6 

59.9 

67.3 

72.2 

72.1 

70-608 

44.9 

47.7 

52.4 

62.7 

69.7 

72.4 

60-508 

33.1 

33.3 

33.5 

39.1 

53.0 

69.4 

50^408 

30.2 

30.4 

30.6 

30.8 

31.0 

51.0 

40-308 

26.5 

26.7 

26.9 

27.2 

27.4 

27.7 

30-208 

27.5 

27.8 

28.0 

28.3 

28.5 

28.7 

20-108 

29.0 

29.2 

29.4 

29.6 

29.8 

30.0 

10-08 

27.4 

27.6 

27.9 

28.1 

28.3 

28.5 

0-10S 

26.1 

26.4 

26.6 

26.9 

27.1 

27.4 

L0-20S 

25.0 

25.3 

25.5 

25.8 

26.1 

26.3 

20-30S 

27.9 

28.1 

28.4 

28.6 

28.8 

29.1 

30-406 

30.0 

30.2 

30.5 

20.7 

30.9 

31.3 

40-50S 

35.2 

35.3 

35.5 

35.7 

38.1 

53.8 

50-60S 

44.0 

44*5 

46.1 

51.3 

60.6 

72.3 

60-70S 

65.5 

65.8 

67.1 

71.0 

72.5 

72.3 

70-80S 

74.2 

74.0 

73.9 

73.8 

73.6 

73.5 

80-90S 

75.4 

75.3 

75.1 

75.0 

74.9 

74.7 


- 48 - 







4PPBTOH B 


Zonal Variations <•><* fitted «nd 
Reflect -orf Rwdi » n t fluxes 
with Particulate Loading 

By operating on the data lirted in Appendix A with equations (1) and 
(2) it was possible to generate the zonal variation in emitted long-wave 
radiation (V e ) and reflected short-wave radiation (W r ) as a function of 
particulate loading. 

The results from these calculations are reported below. 




TABLE Bl: The Zonal Variation of Emitted Radiant Flux 

with Particulate Loading for January (W 0 in vatts/m ) 


Latitude 

Particulate Loading 

Zone 

.300 

•323 

.346 

— .369 

.392 

.as 

90-8GH 

163.0 

161.6 

158.7 

154.2 

147.1 

135.1 

80-70H 

169.8 

167.8 

165.1 

161.2 

154.4 

151.5 

70-60N 

183.4 

181.7 

178.4 

173.9 

165.2 

152.6 

60-5GH 

198.2 

196.7 

194.7 

189.9 

181.8 

166.1 

50-40H 

213.1 

211.9 

210.0 

206.2 

199.9 

185.9 

40-30H 

231.3 

230.3 

228.7 

225.1 

219.3 

207.1 

30-20K 

245.9 

245.0 

243.6 

240.3 

235.0 

223.6 

20-1GH 

252.2 

251.5 

250.0 

246.8 

241-9 

231.7 

IO-OR 

253.3 

252.5 

251.0 

248.1 

243.0 

233.5 

0-10S 

254.5 

253.8 

252.4 

249.4 

244.6 

234.7 

10 -aos 

255.7 

254.7 

253.6 

250.3 

245.7 

236.1 

20-30S 

254.9 

253.9 

253.0 

249.8 

245.1 

235.4 

30-406 

244.3 

243.4 

242.5 

239.1 

234.3 

224.2 

40-50S 

227.8 

227.4 

226.3 

223.1 

218.0 

207.4 

50-60S 

211.7 

211.4 

210.5 

207.0 

205.4 

192.0 

60-7DS 

198.7 

198.6 

197.5 

194.6 

190.3 

182.0 

70-80S 

191.7 

191.4 

190.1 

187.3 

183.0 

174.4 

80-90S 

178.6 

178.4 

177.3 

174.6 

171.2 

163.1 






TABLE B2: The Zonal Variation in Quitted Radiant Flux 

with Particulate Loading for July (W„ in watts/m^) 

V 



Particulate Loading 


.300 

.323 

.346 

.369 

•392 

.415 

90-80N 

192.4 

191.2 

189.3 

184.8 

177.8 

166.5 

80-70N 

201.9 

200.6 

197.9 

192.7 

184.5 

171.8 

70-6CW 

216.1 

214.3 

210.8 

203.7 

194.1 

179.5 

60-50H 

227.6 

226.1 

223.6 

217.6 

207.6 

190.0 

50-40H 

240.5 

239.4 

237.4 

232.6 

225.2 

208.1 

40-30H 

253.4 

252.7 

251.2 

246.8 

240.5 

227.0 

30-20K 

259.2 

258.0 

256.8 

252.9 

247.0 

235.6 

2O-10M 

253.2 

252.5 

251.2 

248.1 

242.7 

232.6 

liBgiMIfegfe 

252.6 

252.0 

250.4 

247.4 

242.4 

232.9 

p-ios 

253.9 

252.1 

252.1 

248.9 

244.0 

234.2 

illtwMcJBW 

253.3 

252.7 

251.7 

248.7 

243.6 

233.9 

20-30S 

245.7 

245.2 

244.0 

241.1 

236.3 

225.8 

30-40S 

232.7 

232.2 

231.0 

227.9 

222.9 

211.9 

40-50S 

214.8 

214.6 

213.3 

209.9 

204.4 

190.7 

50-60S 

199.5 

199.2 

197.8 

193.4 

186.6 

172.0 

60-70S 

184.1 

183.6 

182.2 

177.5 

170.9 

158.7 

70-80S 

171.3 

171.0 

169.5 

165.2 

159.0 

147.4 

80-90S 

156.0 

155.7 

154.6 

150.6 

144.5 

133.5 







TABLE B3: The Zonal Variation in Reflected Radiant Tlux 

with Particulate Loading for January (W r in vatts/m*) 


Latitude 

Particulate Loading 

Zone 

.300 

.323 

.346 

.369 

.392 


90-80N 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

80-7QN 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

70-60N 

12.46 

12.54 

12.71 

13.17 

13.51 

13.49 

60-50N 

36.09 

37.50 

39.72 

44.47 

48.47 

53.95 

50-40N 

46.59 

47*44 

50.29 

56.68 

66.76 

89.35 

40-30N 

62.19 

62.61 

63.04 

63.46 

63.88 

81.44 

30-20N 

77.97 

78.53 

79.37 

79.93 

80.49 

81.32 

20-10N 

87.97 

88.66 

89.68 

90.70 

91.37 

92. 

10-0N 

102.0 

102.8 

103.9 

104.7 

105.9 

106.7 

0-10S 

116.3 

117.1 

118.4 

119.3 

120.2 

121.5 

10-20S 

123.6 

125.0 

125.9 

127.3 

128.2 

129.2 

20-30S 

130.8 

132.3 

133.3 

134.7 

135.7 

136.7 

30-40S 

134.0 

135.4 

136.4 

137.9 

138.9 

139.9 

40-50S 

146*8 

149.8 

150.7 

151.7 

152.7 

202.3 

50-60S 

165.0 

167.4 

174 1 

194.0 

227.8 

299.1 

•S0-70S 

241.2 

243.1 

2 c j.5 

270.0 

300.2 

332.8 

70-80S 

333.4 

333.4 

336.7 

347.6 

347.2 

346.7 

80-90S 

365.6 

365.1 

364.6 

364.2 

363.2 

362.7 






TABLE B 4 : Hie Zonal Variation in Radiant Flux with 

Loading for July (W r in watts/m^) 


Latitude 

Zone 

Particulate Loading 

.300 

.323 

.346 

.369 

.392 

•415 

90-80N 

300.5 

311.4 

330.3 

338.3 

337.8 

337.8 

(80-70N 

248.2 

258.7 

273.8 

307.6 

330.0 

329.5 

70-60N 

200.4 

212.8 

233.8 

279.8 

311.0 

323.1 

60-50N 

150.3 

151.2 

152.1 

177.5 

240.6 

315.0 

50-4QN 

140*4 

141.3 

142.2 

143.1 

144.1 

237.0 

40-30N 

123.6 

124.5 

125.4 

126.8 

127.8 

129.2 

3O-20K 

126.1 

127.5 

128.4 

129.8 

130.7 

131.6 

20-10N 

127.2 

128.0 

128.9 

129.8 

130.7 

131.5 


111.7 

112.5 

113.7 

11.45 

115.4 

116.2 

p-ios 

95.91 

97.01 

97.94 

98.85 

99.59 

100.7 


19.13 

80.08 

80.71 

81.66 

82.61 

83.25 

20-30S 

72.37 

72.89 

73.67 

74.18 

74.71 

75.48 

30-40S 

58.82 

59.22 

59.81 

60.20 

60.59 

61.38 

;o- 5 os 

45.11 

45.24 

45.49 

45.75 

48.83 

68.95 

50-60S 

29.21 

25.55 

30.61 

34.06 

40.23 

48.00 

y0-70S 

10.11 

10.16 

10.36 

10.96 

11.19 

11.16 

70-80S 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

J0-90S 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 






APPENDIX C 

Variations In Global fl"d spheric 
Weighted Averages of Surface Tarnnarature and 
Emitted «nd Reflected Radiant Fluxes 
With Parti ««1 «+• a Tiding 

Using the data from Appendicies A and B, global and hemispheric 
weighted averages for surface temperature, emitted radiation and reflected 
radiation as a function of particulate loading for January and July were 
generated using equations (3) and (4) • 

These data are reported below. 
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TABLE Cl; Global and Heed. spheric Weighted Averages of Surface 

Temperature with Particulate Loading for January (Degrees K) 


Particulate Northern 

Southern 


.ceding 

Hemisphere 

Hemisphere 

Global 

.300 

284.0 

290.5 

287.2 

•323 

283.3 

290.2 

286.7 

.346 

282.4 

289.5 

285.9 

.369 

280.3 

287.8 

284.1 

.392 

276.8 

285.1 

281.0 

.4L5 

269.8 

279.6 

274.7 

TABLE C2: 

Global and Hemispheric Weighted Averages of 
Surface Temperature with P*rticulate Loading 


for July (Degrees K) 


particulate Northern 

Southern 


fading 

Hemisphere 

Hemisphere 

Global 

.300 

293.6 

285.0 

289.3 

.323 

293.0 

284.7 

288.9 

.346 

292.1 

284.0 

288.0 

.369 

289.7 

282.0 

285.9 

.392 

285.9 

279.0 

282.5 

.415 

278.9 

272.3 

275.6 
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TABLE C3: Global and Hemispheric Weighted Averages of 
Quitted Radiant Flux with Particulate Loading 
for January (W Q in w/m 2 ) 


particulate Northern 

Southern 



HemisDhere 

HemisDhere 


MM 

228.0 

237.6 

232.8 

.323 

227.0 

237.0 

232.U 

.346 

225.1 

235.8 

230.5 

.369 

221.6 

232.7 

227.1 

.392 

215.5 

228.3 

221.9 

•435 

204.1 

218.2 

211.1 

TABLE C4‘ 

Global and Hemispheric Weighted Averages of 


Quitted Radiant Flu 

x with Particulate Loadin 


for July (W a In w/m 
© 



particulate Northern 

Southern 

1 

Loading 


HemisDhere 


mm 

244.0 

228.9 

236.4 

.323 

242.9 

228.3 

235.6 

.346 

2a- 1 

227.2 

234.2 

.369 

236.8 

223.7 


.392 

2 * 2.1 

218.3 

22' .2 

.415 

217.4 

207.0 

212.2 
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TABLE C5: Global and Hemispheric Weighted inrerages of 

Reflected Radiant Flux with Particulate 
Loading for January (W r in w/m^) 


particulate Northern 

Southern 

GLobel 

1 

Hemisohere 

Hemisohere 


.300 

63.96 

154.0 

109.0 

.323 

64.63 

155.4 

110.0 

.346 

65.77 

157.5 

111.6 

.369 

67.52 

162.4 

115.0 

.392 

69.66 

168.7 

119.2 

.415 

75.93 

185.0 

130.5 

TABLE C6: 

Global and Hemispheric Weighted Averages of 
Reflected Radiant Flux with Particulate 


Loading for July (W^ 

r in v/m 2 ) 


particulate Northern 

Southern 


joadine 

Hemisohere 

Hemisohere 

Global 

.300 

141.3 

59.02 

100.1 

• 323 

143.7 

59.56 

101.6 

.346 

147.0 

60.15 

103.6 

.369 

.392 

155-4 

165.7 

61.06 

62.50 

108.2 

114.1 

.415 

186.1 

66.27 

126.2 
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APPENDIX D 



The following computer program was used to generate the data 
which appears in Appendicies 6 and C. It also accomplishes the 
conversion of units from kLys to w/m . The program is written in 
standard FORTRAN IV coding. 
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r 


DIMENSION OS( l8,l2).HR«>A( 18*6 >*RNJU( 18.6) »TOJA(l 8.6) » TO JU< 18*6 >» 
1ALRJA( 18*6 )*ALDJU(16*6>/ENJA(t8*A)*ERjU(lA*6)*REFjA( 18*6)* 
IR?TjU( T8*6 ) .QS«(1 8*12) » *N JAN ( 18 * A ) *RN JUN (18.6)# 

1T0JAK( 18»6)*T0JUK( 18*6) 

REA0(5*100)( (OS( I, J)»I«l,ltt), J«1 .i?),((RNJA(I.J)*I»1.18)»J»1*6)* 
1 ( (RNjUI I»J)*1*1*18)*J*1#6)* ( ( T()JA( I»J)*I*t *18)*J*1*6)* 
1((T0JU(1*J)»I*1*18>*J*1#6).((AI BJA(i»J)*I*l*18)*J*l*6)* 

1 ( ( ALB JO ( I » J ) # I s W 18 ) * J*1 *6 ) 

WRTTE(6*99)QS 

NRf Tt( 6* 100)RNJA#RAJU* TUJA* fOJU. Al 8JA.ALBJU 
99 FORMATl 1'2'F 10.3) 

100 format (6Fio.3) 

c Ja aJu stand For January t July respectively* er stands For emitted 
C RADIATION* HEF REFERS 1 □ THt reflected radiation* 

CF>4*1 36E*0 7/ ( 66*00 * *X0 • ) 

DO 10 J*l*12 
DO 10 1*1*18 

10 QSM( 1 • J)*QS( I » J)*CF 
RRlTE (6*101) 

101 Format ^j L T50.' , 0S(NAUS/MtTER**?)^) 

DO IT T*1 * 9 

K. 90 ->(1*10) 

I«K*l0 

11 MRITE(6*10?)L*N*(QSn(I»J)*J*1*12) 

102 FornATC *# 12* "*"• 12» ,, N ,, »12F 10*31 
00 12 T*10» 18 

K«(I-10)*10 

L*K»10 

12 MR|TE(6*103)K*L*(uSm(I*J)*J*1*1?) 

103 format ( * "*i2* ,, *' , *i2."!>"»i2Ho.3i 
Do 13 j*l*6 

00 13 1*1*18 

13 RNuAN(I*j)«knjh(I*u)*CF 

Nrt TE( 6* 10* ) 

io* format ( w , *.T3 o»"rn*jan imatts/meter**?>") 

DO 1* 1*1*9 
K»90"( 1*10) 

L*K*10 

1* HR)TE(6*105)L#R*(RNJAN'11*j)* j«| .a) 

105 Format i * ".i2."*".12. ,, n ,, »6Ho.3> 

DO 15 1*10* 18 

A*(I-10)*10 

L*K»l0 

15 NRTTE(6*106)K*L*(RNJAN(1»J)»j«1 *A1 

106 FORMATl* **;2. M iSl2*"5"*6Fro*3) 

00 16 J*1 *6 

DO 16 1*1*18 

16 RNJUN( I * J)*RNjU( I • J )*Cl 
MR I TE( 6* 107 ) 

107 FORMAT ( w *»T30 *"HN“JULY(kaTTS/METER**2)") 

DO 17 1*1*9 - - 

K*90*( 1*10) 

L*K ♦ 1 0 

17 NRlTE(6*)08)L»A*lRNJUNt l*J)» J*1 .A) 
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106 FoRMATt" 

~ “ Do 18 T»10*18 

K«(I-10)*10 
L.KMO 

16 NRlTE(6*109)K*L,(RNJU*U*J)#J.l.Al 
“109 ~T0RHATr r 

00 31 J*l»i 
00 31 t«l • 1 6 

31 T0JAK(l,J)*T0JA(I»J)«2/i,0 

MRlU(6»130) 

130 _Fj)RHAT(" "»T30**T0“JAl0t(iHEtS K)») 

' 00"3? !»1*9 

K«90*( 1*10) 

l*k»lo‘ 

32 NRtU(6»131)L*N*CT0jAKU*J)»J«l v A) 

T3l ' FoRHAtt" "•T2*"*"*l2* ,, N''*6F10.3) 

On 33 1 * 10*16 

K«Cl*I01*iO‘ 

L*K*10 

" JT ' WIFfTEC 6 #~i 32 )R *L •'( TO jATOTI •J)»J*t»A) 

132 FoRMATt* •*I2»_*- ,, *l2*"i ,, »6F10«31 
Dfi - 34' J*l*6’ 

D0__34 U l»16 

3 « 'ToJUKCI*J)*T 6 JUfi#J> 4 * 27 j.O 
MR1 TE(6»133) 

“ni “ Format! * "*T 3 o»"To*Ju(tfLGREts *>•> 

00 35 I*l»9 
K»90“l 1*10) 

L*K +10 _ _ 

35 RRlTE(. 6 il 34 )L*k#CTUJUKU.J)*J»l • A) 

134 FoRMATC 1 "*I2* ,, - , **1?*' , N"*6F10*3) 
f DO 36 1*10*16 
I K«(I- 10)»10 

j L«K4lO 

! 36 MNfTE(6*135)K*L»( TO JliK( i»J)» v *l • A 1 

f 1*5 TnRMATC *»i2. Hi "»l2»'*W6FlC*3> 

[ on 19 j* i » 6 

, DO lv I«l»lt) 

19 ERJAU*J)*I .97*05*1 J*1>*C1.-ALBJAC_I*J)>)-RNJAN<I»J) 

X OSh( I • 1 ) IS THE LAOITuDINAL VARIATION OF THE SOLAR CONSTANT IN THE 

C HO*fH OF JANUARY . 

r kRlUl6.ilO) 

| no Format t * *.t 20 *"Emitteu radiati. 'Januahy»*> 

• DO 20 1 • 1 * 9 

r - ' K*90-(I*10) 

L*K+10 

I — 20 NRfTE<6*lll)L*F.CERjAU*j),J«1.6> 

in Format c* "* 12 *' , »"» 12 * ,, n , *,61 io.3i 
Do 21 l*lo» id 
K«CI-10)*10 
I L»a*10 

21 NR1TE(6*U?)K*L»UKJAU*J>*J«1*6) _ 

1 12 FDrHATC* •*# 12» Wmm » r2» w S w ,6F 10,3) 

DO 22 J» 1 #A 

DO 22 1*1* 18 

22 ERJU(l*J)*(*97*eSN(I*F)*(l»-ALBJII(I*J)n*RNjUM(|»J) 

C QSN( 1 ,7) IS THE LADllUUINAL VAH 1 AT I UN OF THL SOLAR CONSTANT IN THE 
C MONTH OF JULY 

W*rTEc 5 ill 3 > 

113 Format i * "*t 20 **»ehitteu radiationc jult j") 

DO 23 1*1*9 

K«90*( 1*10) 


JK^fAGE 


OP POOR 


IS- 
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L*K*10 

23 WRITE! 6*1 14)1 ***(EhJUil*J)»J*l *6) 

114 FORMAT!" "*l?."«".l2*"N",6F 10.3) 

DO 24 1*10*16 

*!»( 1*1 J)*10 
L*k*lt> 

24 WR|TE(6*lt5)k.L*(F*'JUU*J)*J« *4) 

115 FORMAT!" ». 12. "S",6fl0. !) 

DO 25 J* 1 » 6 

On 25 1*1*18 

25 REFJA(l*J)»*97*oSwi I » 1 ) *Al rt JA i [ »J> 
NR|TE(6*Vl6> 

116 FORMAT!" "*T2Q-"HE» LtCItO HAD! ATTnN! JANUARY ) 
DO 26 1*1*9 

K«90*( 1*10) 

L*A *10 

26 HRITL(6*U7)L*a*1KLFJA(1,j)*J>1.a> 

117 FORMAT!" "*I2*"-".I2»"N",6F 10*3) 

DO 27 1*10*18 

K*Cl-10)*10 

L*KMO 

27 WRTTL(6*118)H#L # (REF JA! i * j ) • J sl ,*) 

118 FORMAT!" "*12."-".I2*"5",6F10*3) 

DO 28 J*i »6 

DO 28 1*1*18 

26 REF JU(l»J)**97*uSw( I»7 )*aLBJU( I*J) 

WRI TE (6*119) 

119 FORMAT!" "*T20*"REFLE(.Ttn HAOIiTTONtJULT >« 
DO 29 1*1.9 
K«90-( 1*10) 

L*K*10 

29 WRITE! 6»l20)L»K*(REf JUU.JJ.j*! *6.) 

120 FORMAK" 

UO 30 1*10*18 
K*< I“kO)*lO 
L*K ♦ 1 0 

>0 «RlTLit.l21)k*L.lRt* JUti.J)*^* .6.) 

li t Format t " »*i?*"«"*l2**s",6U0. o 

C » INj GIUOAL & HEMISPMEH iCAL WLlUhTu AVERAGES 










0? 


C JANUARY* buRFACt UmP. n JULY. NF 

c jult* emItteo radiation in January 

C IN JANUARY 4 JULY 
MR!Ttl6*122) 

122 FORMAT!" "*T30."AVC«A!»L 

Call aver(Tojak) 

WRITE (6.123) 

123 FORMAT!" "*T25* "AVtwAOL 

call averctojuk) 

WRITE (6*124) 

124 FORMAT!" ", T25. "Avt < AUL 
CALL AVLR(RNJAW) 

WRI TL ! 6 * 1 25) 

125 FORMAT!" "*T25* "AVL YAUt 


' RADIATION 
ft JULY AND 


for Surface temp. in 

IN JANUARY & 
REFLECTED RAOIATlUN 


SUKFACc TEMPERATURE FOR JANcAHY") 


Sure ac * tfmperaturE for July") 


net rahatIun Fur januart*) 


net rajiatiun for July") 


Call aver(knjiin) 

WR I TE ( 6 » 1 26) 

126 FORMAT!" ",T25»"AVt <Al»t EMlTTFl) RADIATION FOR JANUARY") 

call avercerja) 

WRlTt(h»127) 

127 FORMAT!" ",T25» "AVE yAGt EmITTEu RADIATION FOR JuLY") 

Call avlr(Ehju) 

WRI I t ( 6* 1 28 ) 

126 FORMAT!" "*T25*"AVtRAUt REFLECTED RADIATION FqR JANUARY") 
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Call aver!REFja> 

WRTTL(6',129) 

129 F 0 rHAT!" "*T25*"AV£RAGL RETi-ECTEO RADIATION FOR JULY") 
Call avEr!REFju> 

STOP 

End 


Subroutine aver (m 

REAL north 

01 MENS! ON X( 18*6)»SUHNtO)#NUHTM(A)» SUMS! 6 ) • SOUThC 6 ) .GLOBAL! 6 ) 
On 16 J*1 »6 
*ETGHT»0.0 
H«0 .0 " 

SUHNIJ)«0.0 

Oo' iR 1 * 1 #9 

NEIGHT*»IEIGHT»!SIN!3.1*/2.-H>-SIN(C8.*3.1A/18. )-H)) 
"5UMN(J)«SUHN( JTt'(X(i»J>*( SIN( 3.1 A/2."H)“SlN((d**3*U/l»* )•«))) 

19 H«H*!3»14/16»> 

U0RTH('j)1»5uHNt J)/NEltiHT ' 

NRtTE(6»101)N0RTH(J) 

1 01 ForIATC" "*"ThE NORTHERN HEMISPHERICAL WEIGHTED average for 
1 THE QUANITY NOTED ABOVE 1S".F10.3> 

is Continue 

on 20 J»1»A _ . 

RE16HT»0i0 

HsO.O 

Suhsi J)*0.0 

00 21 1«10> 16 

»lFTUHT»NEICHT*tSlNt 3.1^/18.+H)-STN(H) ) 

SuhSI J)«SUMS( J)*(X! I» J>*(SIN( 3. 1 A/18.*H)-SlN(H) ) > 

'21 H*h*1 3*14/ 18.) 

SouTHC J)*SLHS( J) /HEIGHT 
NrtTEC6»10?)S0UTH( j) 

102 FORMAT!" "."The SOUTHERN HEmISPhFRICAL mETGhTEO average for 
1 The 0UANITY NOTED ABOVE IS H .F10.3> 

20 continue 

DTT 22 J* 1'. A “ 

GLOBAL! J)b(SUHN( J) ♦SUMS l j ))/(?. *wF I GHT ) 

NR|TE(6#103)GLOBAL( J) 

103 FORMAT!" "."THE GlUBAL "EIGHTFO AVERAGE FnR ThE GUAnITY NOTED 
1 AB0VE"*F10.3) 

22 CONTINUE 
RETURN 
end 
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